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Abstract 

The spatial distribution of interactions in highi energy collisions of heavy nu- 
clei is discussed using the wounded nucleon, binary collision, hard sphere, and 
colliding disk parameterizations of interaction densities. The mean radius, its dis- 
persion, and the eccentricity of the interaction region are calculated as a function 
of impact parameter. The eccentricity is of special interest for comparison to mea- 
surements of anisotropic flow. The number of participants and binary collisions is 
also tabulated as a function of impact parameter. 

1 Introduction 

In this note we discuss the spatial distribution of initial interactions in the high en- 
ergy collisions of heavy nuclei, in order to gain some insight into geometrical effects 
in such collisions. We apply a simple, widely used formalism (e.g. ^) that is 
often attributed to Glauber to calculate the density of interactions projected onto 
the plane transverse to the beam axis for four commonly used weighting functions: 
number of wounded nucleons, number of binary collisions, colliding hard spheres, and 
colliding disks. The wounded nucleon and binary collision calculations incorporate 
the usual Woods-Saxon parametrization of nuclear density. We compute the mean and 
dispersion of the radius and the eccentricity of the interaction zone as a function of 
impact parameter for symmetric collisions of heavy nuclei. Using the same calcula- 
tion we tabulate the number of participants for the wounded nucleon and hard sphere 
paramterizations, together with the number of binary collisions, as a function of impact 
parameter. 

The distributions calculated are implicit in all model calculations of low or high so- 
phistication that incorporate similar parameterizations of nuclear geometry, and in that 
sense nothing new is presented here. However, we have not found these distributions 
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collected and discussed explicitly in the literature, and that is the the purpose of this 
note. A more detailed discussion of these and related results was presented previously 




The formalism is described briefly in section 2. Section 3 discusses the transverse 
profile of the interaction volume as a function of impact parameter. Section 4 presents 
a parametrization of the eccentricity of the interaction region for non-central collisions. 
Section 5 tabulates the number of participants for the wounded nucleon and hard sphere 
paramterizations, together with the number of binary collisons, as a function of impact 
parameter. 



2 Formalism 

The coordinate system is defined in Figure |l]. We utilize the standard nuclear thickness 
function T A 

Ta(|s1)= J dzpA{z,s). (1) 
For the nuclear density we use a Woods-Saxon distribution, 

PA{r)^Po-- / p, - , , (2) 

1 -I- e(''~-^A)/a 



where r = Vs^ + z^, Ra = 1.12 • k^'^, and po = 0.159 GeV/W and a=0.535 fm 
for ^^^Au. PA is normalized so that J (i^sTA(|s|) = A. 

In Fig. |l|, the nuclear impact parameter is denoted b, and the distance to a point in 
the transverse plane from the center of either nucleus is denoted bA and be- The vector 
from the origin to this point is written 

s = bA - = be + ^i. (3) 

where x is the unit vector in the x direction. 

To calculate the density of interactions as a function of impact parameter for a given 
process, we utilize the limiting cases of high and low nucleon-nucleon cross section for 
that process: 

• Wounded Nucleon Scaling: We consider the number of nucleons at s that are 
struck at least once by the nucleons in the oncoming nucleus, where "struck" 
means inelastically excited with nucleon-nucleon collision cross section||] ctnn- 
In the transverse projection, the density of wounded nucleons per nuclear colli- 
sion is given in units 1 / fm^ by: 



= TA(bA)-(l-e-TB(bB)<TNN) + TB(bB)-(l-e-TA(bA)CTNN)^ (4) 



ds2 

where the dependence on nuclear impact parameter is via Eq. |[ 



'(TNN~ 30 mb at the SPS (y/s^=2a GeV) and 40 mb at RHIC (v^s5JS"=200 GeV). 
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Figure 1 : Transverse coordinate system. A and B are the masses of the colliding nuclei. 
The incoming beams are (anti-)parallel to the z axis. The x-z plane is commonly called 
the reaction plane. 



• Binary Collision Scaling: More generally, we mean those processes with suffi- 
ciently small nucleon-nucleon cross section ctJ^n'^ '^^at (Ta-ctJ^n^)^! ^PPly 
the label "hard" because we are usually refering to high momentum transfer pro- 
cesses, or hard scattering). Nucleon-nucleon interaction probabilities can then be 
summed for the total interaction probability, and the number of hard scatterings 
per nuclear collision goes as the number of binary nucleon-nucleon collisions. 
The density in the transverse projection of hard processes per nuclear collision, 
in units 1/fm^, is then 

^^^ = fTh-d.T^(bA)-TB(bB). (5) 

Integrating over the transverse plane (and changing integration variables), the 
total number of hard processes per nuclear collision is given as a function of 
nuclear impact parameter b by[[l|, ^: 

Nhard(;^) = ah'^rJ'^ • / dhTA{\s\)TB{\b- s1) = a^-d • Tab(&) (6) 



3 Transverse Profile of the Interaction Volume for Non- 
central Collisions 

The azimuthal anisotropy of momentum distributions in nuclear collisions can be re- 
lated to the orientation of the event plane (the x-z plane in Figure [T]) for noncentral 
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collisions over a wide range of energies [|5[|, and has been quantitatively studied up 
to the highest energy nuclear collisions at the SPS |^ and RHIC The lowest or- 
der harmonics are referred to as directed and elliptic flow, and are characterized by 
coefficients vi and V2 respectively in a Fourier expansion of the azimuthal angle or 
momentum distributions . 

Elliptic flow at midrapidity has received particular attention recently because of a 
possible connection to the Equation of State (see discussion and references in [^). To 
help distinguish dynamics from purely geometrical effects, it has been suggested [jsj ^ 
that the measured V2, the elliptic anisotropy, be scaled by the eccentricity of the reaction 
volume. This is defined to be [|[ ||] 

^ <y'^> ~ <x^ > _ _b_ 

^ ~ < Z/2 > + < a;2 > 2Ra 

where < . . . > indicates the spatial average over the transverse plane weighted by a 
density such as that of wounded nucleons (equation^. The approximation is the ratio 
of the lengths of the axes of the overlap region in Figure [l|, {y \x=o)/{x |y=o), not 
weighted by nuclear density [|| . 

We calculate the transverse density distribution of the reaction volume as a function 
of impact parameter, utilizing four different weighting functions: 

• Wounded Nucleons: The transverse density profile is calculated using a Woods- 
Saxon density distribution and Eq. Q appropriate for the bulk of particle produc- 
tion. We use (Tnn=30 mb, but also investigate the sensitivity of the computed 
quantities to this parameter 

• Binary Collisions: The transverse density profile is calculated using a Woods- 
Saxon density distribution and Eq. ^, appropriate for hard processes such as jet 
and J /'J' production. 

• Hard Sphere: The transverse density profile is calculated for colliding sharp- 
edged spheres, with the density defined as Ta+Tb. This corresponds to the limit 
of the Wounded Nucleon density in which the Woods-Saxon parameter a in Eq. ^ 
is small and (Jnn in Eq. ^is large. The radius of the hard sphere corresponding 
to an Au nucleus is 7.24 fm, increased from the Woods-Saxon value of 6.52 fm 
in order to obtain the same total interaction cross section. 

• Two Dimensional: The transverse density profile is simply the area of overlap 
region in Fig. [l|, with uniform weighting (i.e. without taking into account the 
nuclear thickness in the z-direction and corresponding variation in density when 
projected onto the plane). The radius corresponding to a Au nucleus is 7.24 fm, 
as in the Hard Sphere calculation. 

We first present the transverse density profiles graphically, and then calculate mo- 
ments to compare the distributions quantitatively. Figures || to ^ show the transverse 
density profiles for AuH-Au collisions as linear contour plots in one quandrant of the co- 
ordinate system defined in Fig. |l|, for the four weighting functions at different impact 
parameters. The Wounded Nucleon and Binary Collision profiles have very similar 
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Figure 2: Transverse density of Wounded Nucleons as a function of (x,y) (see Fig |l]), 
for collisions of Au+Au at impact parameters b=0, 4, 8 and 12 fm. 
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Figure 3: Transverse density of Binary Collisions as a function of (x,y) (see Fig gj), for 
collisions of Au+Au at impact parameters b=0, 4, 8 and 12 fm. 
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Figure 4: Transverse density profile of interacting hard spheres as a function of (x,y) 
(see Fig |l]), corresponding to Au+Au at impact parameters b=0, 4, 8 and 12 fm. 
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Figure 5: Two dimensional transverse densityprofile as a function of (x,y) (see Fig[l|) 
corresponding to the region of overlap in Fig |l| for Au+Au collisions at impact param- 
eters b=0, 4, 8 and 12 fm. 
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Figure 6: Mean and dispersion of radius of the interaction region as a function of 
impact parameter for Au+Au collisions, for the four weighting functions. The units of 
both axes are fm. 



shape, with the Binary Collision distribution falling off slightly faster from the ori- 
gin. The less realistic Hard Sphere and Two Dimensional calculations exhibit larger 
aspect ratios for noncentral collisions and unphysical sharp declines in the density at 
the boundaries of the overlap region. 

Figure ^ shows th e mean transverse radius (< |s| > for Fig. [l|) and its dispersion 

(y/<|s|^> -<|sl>^) for the four weighting functions, as a function of impact pa- 
rameter for AuH-Au collisions. As is seen in the density profiles themselves, the Two 
Dimensional and Hard Sphere functions give larger mean radii than the Wounded Nu- 
cleon and Binary ColUsion functions. The dispersion in the radius is similar for all 
four functions, and is a weak function of impact parameter The values at large impact 
parameter are dominated by the treatment of the nuclear surface. 

Note that in the case of Binary Collision weighting, the mean (and median) radius 
is about 3.5 fm for the most central collisions. In other words, about one half of all 
produced e.g. J/^ or jets are generated farther than 3.5 fm. in the transverse plane 
from the center of the reaction zone in central collisions, with a distribution in radius 
which has a half- width of about 1 .5 fm. Rather few of the J/^ or jets are produced near 
the center of the reaction, due simply to the geometry of nucleus-nucleus collisions. 

Figure ^ shows the eccentricity (Eq. ^ of the interaction region as a function of 
impact parameter for AuH-Au collisions, for the four weighting functions. Also shown 
is the approximation b/(2RA) from Eq. M. The eccentricity of the Wounded Nucleon 
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Figure 7: Eccentricity (Eq. |7|) of the interaction region as a function of impact parame- 
ter for Au+Au collisions, for the four weighting functions. 

and Binary Collision models are similar and significantly smaller than those of the 
other models. Parametrization of e in the Wounded Nucleon Model and its dependence 
upon ctnn is given in the next section. 

4 Parametrization of e and S in the Wounded Nucleon 
Model 

In this section we give a more detailed discussion of the characterization of the geom- 
etry of the interaction region in the Wounded Nucleon Model. The elliptic anisotropy 
of the overlap region e, defined in Eq. 0, is shown as a function of impact parameter 
in Fig. 0for all four models for AuH-Au collisions, e is parametrized for the Wounded 
Nucleon Model ((7nn=30 mb) as 

e(b) = -0.0469a; + 2.754x2 - 4.797a;^ + 4.852a;'' - 2.492a;^ (8) 

where x = b/(2R) and 2R 15 fm. for Pb-nPb collisions and 2R = 14.7 fm. for 
AuH-Au collisions. 

Increase of (Jnn to 40 mb decreases e by 13% at b=2.5 fm and 4% at b=10.5 fm, 
with the change monotonic in b. The decrease of e with increasing ctnn can be un- 
derstood by noting that the value of e is dominated by surface effects. Increasing ctnn 
effectively extends the surface to larger radii, with the consequence that the shape of 
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the interaction region in the transverse projection is less eccentric than for (7nn=30 mb 
at the same impact parameter 

For reference, we also include a parametrization in the same model of the impact 
parameter dependence of the quantity S = nRxRy Here, = ^/<x^>, Ry = 
yj <y'^>, and < . . . > denotes the spatial average used in Equation 0. In [^, S* is used 
to calculate the particle density in the overlap region. The parametrization is 

S'(b)/(7rR2) = 0.164 + 0.0141a; - QMAx^ + 1.026x3 - 0.763a;^ + 0.284x^ (9) 

where, again, x = b/(2R) and 2R = 15 fm. for Pb+Pb colHsions and 2R = 14.7 fm. 
for Au+Au collisions. 



b(fm) 


# part 
(WN, crNN=30 mb) 


#part 
(WN, crNN=40 mb) 


#part(HS) 


#BC (phmiiC^) 


0.5 


368 


376 


389 


2.8 • 10-^ 


1.5 


355 


364 


366 


2.7-10-^ 


2.5 


331 


341 


335 


2.5 • 10-^ 


3.5 


298 


310 


302 


2.2 • 10-^ 


4.5 


262 


274 


263 


1.8 • 10-^ 


5.5 


223 


234 


226 


1.5 • 10-^ 


6.5 


183 


193 


187 


1.2 • 10-"^ 


7.5 


144 


155 


148 


8.5 • 10-3 


8.5 


108 


117 


115 


5.8 • 10-^ 


9.5 


77 


84 


82 


3.7-10-3 


10.5 


49 


55 


54 


2.1 • 10-3 


11.5 


28 


32 


32 


1.0- 10-3 


12.5 


13 


16 


15 


4.0- 10--* 


13.5 


5 


6 


4 


1.3- 10--* 



Table 1 : Impact parameter dependence for Au+Au collisions of the number per nuclear 
collision of nucleon participants in the Wounded Nucleon model (columns 2 and 3) and 
colliding Hard Spheres (column 4), and the number of binary scatterings with cross 
section 1 /xbarn (column 5). 



5 Number of Participants and Binary Collisions 

We conclude with a tabulation of the number of participants and binary collisions as 
a function of impact parameter. The distributions in Figures § to ^ can be integrated 
to calculate the total number of nucleon participants or binary collisions per nuclear 
collision. These are presented in Tables |l] and ^ for Au+Au and Pb+Pb collisions. 
Columns 2 and 3 give the number of Wounded Nucleons (Equation ^ with (Tnn=30 
and 40 mb. Column 4 gives the fraction of volume of the Hard Spheres that interact, 
normalized to the total number of incoming nucleons. Column 5 gives the number of 
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b(fm) 


# part 
(WN, crNN=30 mb) 


#part 
(WN, crNN=40 mb) 


#part(HS) 


#BC inharn-') 


0.5 


388 


397 


410 


3.0 • 10-^ 


1.5 


375 


385 


388 


2.9 • 10-^ 


2.5 


351 


362 


358 


2.6 • 10-^ 


3.5 


318 


330 


322 


2.3 • 10"^ 


4.5 


281 


293 


286 


2.0 • 10-^ 


5.5 


241 


252 


246 


1.6 • 10-^ 


6.5 


199 


211 


206 


1.3 • 10-^ 


7.5 


159 


170 


169 


9.5 • 10-^ 


8.5 


121 


131 


132 


6.7-10-^ 


9.5 


87 


95 


99 


4.3 • 10-^ 


10.5 


58 


64 


69 


2.5 • 10"^ 


11.5 


34 


39 


44 


1.3-10-^ 


12.5 


17 


20 


23 


5.5 • 10--* 


13.5 


7 


9 


8 


1.9- 10-* 



Table 2: Impact parameter dependence for Pb+Pb collisions of the number per nuclear 
collision of nucleon participants in the Wounded Nucleon model (columns 2 and 3) 
and colliding Hard Spheres (column 4), and the number of binary scatterings with 
cross section 1 /ibarn (column 5). 



binary collisions for an interaction cross section of 1 /ibarn (the rate for other small 
cross sections is obtained by linear scaling of this number). 

Tables |l| and || show only a weak dependence of the number of participants on 
tTNN, and good agreement between the number of participants calculated with the Hard 
Sphere and Wounded Nucleon models. The latter point is in sharp contrast to the 
strong difference seen between the models for e (Figure |^), and is due to the fact that 
the number of participants is dominated by the bulk volume, whereas e is dominated 
by the surface overlap. 
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